Staggered chiral random matrix theory 



James C. Osborn 

Leadership Computing Facility, Argonne National Laboratory, Argonne, LL 604-39, USA 

We present a random matrix theory (RMT) for the staggered lattice QCD Dirac operator. The 
staggered RMT is equivalent to the zero-momentum limit of the staggered chiral Lagrangian and 
includes all taste breaking terms at their leading order. This is an extension of previous work which 
only included some of the taste breaking terms. We will also present some results for the taste 
breaking contributions to the partition function and the Dirac eigenvalues. 



I. INTRODUCTION 



II. STAGGERED FERMIONS 
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Staggered fermions are one of the commonly used ways 
to simulate quarks on a lattice due to their relatively low 
computational expense. However they are considerably 
more difficult to deal with theoretically due to the pres- 
ence of extra modes. A single staggered Dirac matrix 
yields four flavors of quarks in the continuum limit due 
to the fermion doubling problem. These flavors are mixed 
at finite lattice spacing and are conventionally referred to 
as "tastes" to distinguish them from regular quark fla- 
vors. The breaking of the SU(4) taste symmetry can be 
reduced by using improved actions, but the taste break- 
ing (TB) is still observed in current simulations and must 
be accounted for when extracting results. 



The form of the taste breaking has been worked out as 
corrections to the chiral Lagrangian by Lee and Sharpe 
with extension to multiple flavors by Aubin and 
Bernard The staggered chiral Lagrangian includes 
all terms of 0{a?) [a is the lattice spacing) that are con- 
sistent with the symmetries of staggered fermions. 



Here we will construct a complete random matrix the- 
ory for the staggered Dirac matrix that incorporates all 
terms of 0{a?). This is an extension of the work in 
Q where only the case of zero topological charge was 
considered, and not all of the terms found in the chi- 
ral Lagrangian could be reproduced in the RMT. This 
RMT can be directly related to the staggered chiral La- 
grangian in the zero momentum limit. We then expect 
to be able to use this model to study the effects of TB 
on low energy quantities such as the partition function, 
chiral condensate, and Dirac eigenvalues. Additionally, 
since the standard method to deal with the extra quark 
modes is by taking the fourth root (or square root) of the 
fermion determinant, we expect to be able to study the 
interactions of the TB with the "rooting" procedure. In 
this work we will study the properties of the RMT itself 
and save the comparisons to direct simulations of lattice 
QCD for later. 



On a 4d lattice the unimproved staggered fermion ac- 
tion, -|~ am]'0, with quark mass m is given by 

1 - 

where C/^ is a set of SU(iVc) matrices representing the 
gauge field and rj^fi = {—1)^"<>^^". Here, for conve- 
nience, we will only consider the case of A^c = 3. The 
RMT given below will also apply for A"c > 3, however for 
Nc — 2 the low eigenmodes are known to be in a different 
universality class, described by the Gaussian Symplectic 
Ensemble [4|. One should be able to extend this model 
to include Nc = 2 in a similar fashion, but we will not 
pursue that here. 

The massless staggered fermion matrix is anti- 
Hermitian and thus has purely imaginary eigenvalues. It 
also possesses a particular form of chiral symmetry 



{T,,Ds} = 0, Fs = (-l)2^'--i 



(2) 



which causes the eigenvalues to come in positive and neg- 
ative pairs: ±iX. It is well known that the above action 
contains doubler modes so that a single staggered fermion 
matrix actually describes four flavors of fermions in the 
continuum limit. An explicit identification of the con- 
tinuum fermions was given by Kluberg-Stern, et al. 
This was done by transforming the 3 color degrees of 
freedom on the 16 sites of each 2^ hypercube into a basis 
of four Dirac fermions (12 components each) on a lattice 
of half the size in each direction. This transformation is 
not unique. In their basis the expansion of the staggered 
fermion operator starts as 



(2a)4 {(7^ (g) l4)Df, + m(l4 I4) - 0(75 (g ^^,5)0 



•}• 

(3) 



The notation {S(E)T) is used to denote the outer product 
of a 4 X 4 spin matrix S and a 4 x 4 taste matrix T with 
==7^! Cm5 = Cm^s, and I4 is a 4 x 4 identity matrix. 
The first part is just the standard Dirac operator for 
four identical flavors with mass m. The remaining term 
is suppressed by a factor of the lattice spacing and breaks 
the SU(4) taste symmetry. There is also a term of 0{ag), 
and the remaining corrections are at least O(a^). 
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In the full theory one usually describes the low energy 
behavior in terms of an effective chiral Lagrangian. For 
staggered fermions this has been worked out to order a? 
and is given by [l|, Q 

C='^{df.Ud^U^)-]-i:om{U + U^)+a''V (4) 

where F and So are the low energy constants related to 
the pion decay constant (with the convention that the 
physical value for F w 131 MeV) and the (absolute value 
of the) chiral condensate respectively. Here and every- 
where below, {X) will stand for the trace of X. There is 
also a mass term for the taste singlet pion (analogue of 
the rj') that we have dropped. 

The taste breaking terms can be divided into two parts 
V = V\t + V2t- The first part contains the single-trace 
terms 

+ C4lY.^^^^^^^5,U) + h.c] 

+ CeY. (^^.U^.^U^) (5) 

and the second part has the two-trace terms 
-V2t=C2yi5][(C,,{/) {^^.U)+h.c.] 

+ ^5A^E[(^A.5t/)(e5Mf^^)]. (6) 

Note that in the original Lee-Sharpe Lagrangian, the two- 
trace terms were Fierz transformed into one-trace terms 
of the form 

- V25 = C2^ [{U^} - {i5U^5U) + h.c] 

+ E [(C^f^^M^^) - {^^M,U^)] ■ (7) 

This is valid in the one flavor case, but not when extend- 
ing to multiple flavors 0]. We will see below that the 
staggered RMT naturally leads to the two-trace form. 



where W is a (A'' + v) x N complex matrix with i/ the 
absolute value of the topological charge and the Dirac 
operator is represented by [1| 

It has been shown that the partition function is universal 
for a large class of weights 0, [1] , but for convenience here 
we take the simplest form of a Gaussian, 

p{W) = exp(-a7V (W^W)) (10) 

with y/a = J^oV/2N {V is the four volume). This model 
is the chiral extension of the Gaussian Unitary Ensemble 
(GUE). For a review of chiral random matrix models, see 
Q and its references. 

By now it is well established that the chiral RMT 
(including all spectral properties) is equivalent to the 
zero-momentum sector of the chiral effective theory [Tol - 
[12} . The equivalence is established through the partially 
quenched partition functions 

= • (11) 

Here Z^^'^"^ is the RMT partition function in the micro- 
scopic limit, defined by taking the limit A^, — > 00 while 
keeping rh = mVl^Q fixed. The chiral effective theory at 
zero momentum (also called the e-regime ^13j ) contains 
just a mass term 

ZlIXa^AA^b}) ^ J dUdet{Ureii''^^+^'^){l2) 

with M — diag({m/}, {mb}). The above expression is 
a supersymmetric generalization of the usual fermionic 
chiral Lagrangian so the determinant and trace must be 
taken to their supersymmetric equivalents, and the in- 
tegration is now over a supersymmetric manifold that is 
compact in the fermionic sector but noncompact in the 
bosonic sector ll]. 

Since the TB terms contain no derivatives, they will 
also contribute to the zero-momentum partition func- 
tion by multiplying the integrand by the extra factor 
exp(— a^yV). Below we will establish the equivalence be- 
tween the partition functions including TB only for the 
simpler fermionic case. In principle, it is necessary to 
show this for the partially quenched partition functions 
as well, in order to establish an equivalence of valence 
quantities including the Dirac eigenvalues. For now we 
will assume that this equivalence holds and that it could 
be obtained from an extension of the proof for the parti- 
tion functions without TB. 



III. CHIRAL RANDOM MATRIX THEORY 

The standard partially quenched chiral random matrix 
theory can be written as 



_ ,„,, ny^idet(Po 

dWp{W) 

nr=i det(Po + mb) 



m/j 



(8) 



IV. STAGGERED CHIRAL RANDOM MATRIX 
THEORY 

To extend the RMT to include taste breaking, we add 
a term proportional to a to a taste diagonal Dirac matrix 



aT 



(13) 
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TABLE I. Mappings from corrections to the chiral RMT [T in Eq. (|13p ] to corrections to the chiral Lagrangian [V in Eq. ([4])] . 
The Gaussian weight is given by exp(— S'r)- The first two types of taste breaking terms are similar to the ones appearing in 
[^. The last terms are new and will generate the two-trace terms. 



where T incorporates the taste breaking terms considered 
below. 

In Q we considered only the case of — 0; fur- 
thermore, the additional terms could only reproduce the 
single-trace terms. Here we will consider the extension 
to 7^ and will also include the two-trace terms. Sim- 
ilar work has been done for the Wilson Dirac operator 

MM- 

We start with the dominant term which is typically 
found to be the C4 term 1,]. For arbitrary v we can 
write it as 

r-E{^^l)^(.. (14) 

where and i?^ are Hermitian matrices of size (N + 
v) X {N + v) and N x N, respectively. Note that this also 
has a chiral and taste structure similar to the leading 
taste breaking term in the expansion in ([31). If we choose 
a Gaussian weight function for these matrices of the form 

eM-mJ2KK) + (Bl)]) (15) 

then one can show that the chiral Lagrangian will get a 
correction term (see Appendix [K\ for details), 

- ^ E (^M5C/CM5f/ + ^,5U^C,5U^) ■ (16) 

Upon equating this to the C4 term in the effective La- 
grangian [and noting that there is a ^5^^ in ([5])], we get 
/3 = aN/2VC4. Note that we require that /3 > for con- 
vergence of the integrals. We could have obtained the 
opposite sign in (1161) if we multiplied (jl4p by i; however, 
this would make that term Hermitian. Thus the sign of 
C4 is determined by the need to have an anti-Hermitian 
Dirac operator in We will discuss this issue more 

in the context of the two-trace terms. 

The C3 term can be handled in a manner similar to 
C4. The Ci and Cg terms can be obtained from matrices 
of the form 



where X is a (N + iy) x N complex matrix. The correction 
to the chiral Lagrangian for this term is given in Table ID 

It was pointed out in ^] that the terms in ([7]) could 
be obtained from a RMT using terms similar to the ones 
above, but the corresponding RMT would contain Her- 
mitian pieces, instead of being strictly anti-Hermitian as 
is the case of the staggered Dirac matrix. By writing 
those terms in the two-trace form one can now find 
a way to add them to the RMT while preserving the 
anti-Hermiticity. 

To do this we need to make linear combinations of the 
terms. For example, we can write the C2V and Csy terms 
as 

^T.(UU + U^)f + ^j:(^^(U-U'))" (18) 
with 

C± = {C2V ± C5v)/2 . (19) 

We can then linearize each of these terms using a 
Hubbard-Stratonovich transformation, such as 

e4(«.(c^+c^^-))' =y'd^e-^[-'"2-(«-(^+^'))] (20) 

where cr is a single real variable and s = \J Cy / 1 Cy \ . 
This term now takes the form of a mass term that mixes 
the tastes (with a mass matrix ^^|Cy jcrs/So). The Cy 
term likewise gives a 75 mass term. The mappings for 
these terms from the RMT to the chiral Lagrangian are 
given in the last two rows of Table ID 

Note again that we now must have Cy < in order 
for this term to be anti-Hermitian. The same condition 
holds for . The coefhcients y are proportional to 
the "hairpin" coefficients 8'^ y which appear in one-loop 
results of chiral perturbation theory 0. The other com- 
binations C\ y do not appear in one-loop expressions 
and therefore have not yet been determined from lattice 
simulations. However, the negative sign for y is con- 
sistent with lattice measurements [iGi] , as are the positive 
signs of all the single-trace coefficients. 
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The inclusion of the two-trace terms in their current 
form may seem a bit ad hoc since they aren't fuh matrices 
hke the other terms, but we can rewrite them in a way 
that seems more natural. As an example, we consider a 
RMT with taste breaking terms of the form 



A 



(fe + c) 






S + (6 - c) (g) Iat 



r (21) 



(A, B are Hcrmitian matrices and 6, c are real scalars) 
with weight 

exp {~(5N[{A'^) + (B^)] - -fNb^ - SNc^) . (22) 

If we make the substitution A' — A + b+c, B' = B + b — c, 
then b and c no longer appear as part of the matrix, but 
appear only in the weight. We can then integrate them 
out to obtain a new weight function. 



exp 



(^-(3N {{[A' - Af) + {[B' - Bf)} 



(23) 



N 



[7(A' + B' f + SiA' - B')^] [1 + 0{1/N)] 



with A' = {A') /[N + v) and B' = {B') /N. In this 
way, we see that the two-trace terms in the chiral La- 
grangian are generated by two-trace terms in the RMT 
potential. One could then consider adding other terms 
such as higher powers of the matrices in the potential to 
reproduce higher order terms in the chiral Lagrangian, 
though we will not pursue that here. 

Ah the terms of the full staggered RMT (SRMT) are 
written out in Appendix [BJ Now that we have the full 
form of the SRMT, we can examine its structure more 
closely. For this, it is convenient to switch to a ba- 
sis where the remnant of chiral symmetry for staggered 
fermions ([2]) is transformed according to 







-In 



'^5 











(24) 



There are several possible choices of basis, all of which 
give a staggered RMT Dirac operator (at m = 0) of the 
form 



(25) 



where i? is a (4A^ -I- 2i^) x (4A^ -I- 2v) matrix. Arbitrarily 
picking one basis gives an R of the form (using the terms 
from Appendix [B]) 




\iA 



3/j. 



-f- d+\(7^ iW + iXi + iXef^„(7f^al 



HB: 



B 



4m 



(26) 



with = (1, -ia*) and d± = ibvp. - CAf^ ± (icy^ ^ ^Ap)- 
Note that since i? is a square matrix, in general for 
nonzero lattice spacing there are no exact zero eigenval- 
ues of the RMT. This agrees with the well-known proper- 
ties of the lattice theory. Also, one can imagine that if the 



taste breaking terms are large enough, then the detailed 
structure of R may not matter, and the low eigenvalues 
are described well by a standard chiral RMT at i/ = 0, in 
agreement with numerical studies (T7l - [20| . Below we will 
take the limit of large taste breaking and show that this 
is indeed the case. 



V. SCALES 

For standard staggered chiral perturbation theory (in 
the p-regime) the size of the taste breaking can be mea- 
sured by the parameter [iGj 



(p) 



Xa2 



87r2F2 



(27) 



where a^A is a "typical" taste breaking term. Taking 
this to be the average pion splitting gives 

A = -^(Ap-f 4Ay-f 6At + 4A^-|-A5) (28) 
16 

where the Ax parametrize the mass shift of the pions 
above the Goldstone pion (the taste pseudoscalar) 



2 2 I 2 A 

m^x = m^p + a Ax 



(29) 



This scale determines the convergence of the taste break- 
ing parts in the perturbative expansion. 

For the zero-momentum chiral Lagrangian (e-regime) 
considered here, the relevant parameter is 



X. 



a'^VC 



(30) 



where C is another measure of the strength of the taste 
breaking, which we will take to be 



C = Ci+AC3 + 4C4 + 6Ce 
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-A 



(31) 



We have ignored the contributions from two-trace terms 
in this definition for simplicity, though one could include 
them if needed. If the scale in ([50)1 is small, then one can 
calculate quantities from a perturbative expansion of the 
zero-momentum effective theory (starting from either the 
SRMT or the chiral Lagrangian) in the taste breaking. 
In this case we expect the low eigenvalues to be nearly 
fourfold degenerate and form clear "quartets." This has 
been seen in lattice simulations with improved actions 

[HHil. 

In the opposite limit, xl^-^ ^ 1, the approximate four- 
fold degeneracy is strongly broken, and the low eigenvalue 
spectrum will resemble that of a single flavor due to the 
remaining unbroken U{1) Goldstone symmetry of stag- 
gered fermions. We will refer to this limit as strong taste 
breaking, and the opposite limit as weak taste breaking, 
independent of the value of ■ 

In typical lattice simulations, the volume is chosen such 
that the lightest dynamical mass stays in the p-regime, 
given by the condition m^ri S> 1 (niTrL « 4 is the usual 
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rule of thumb). In such simulations the smallest eigen- 
values are still typically described by e-regime calcula- 
tions since they can be related to observables with va- 
lence quark masses equal to the eigenvalues. The scale 
below which eigenvalues can be described by the zero- 
momentum Lagrangian (known as the Thouless energy) 
in QCD is given by [M ill 



(32) 



[up to a constant factor of 0(1)]. It is possible for 
eigenvalues below this scale to be in the strong TB 



regime (x^2^ ^ 1), while observables at higher scales, 
e.g. around the dynamical quark mass, exhibit weak TB 
(xl? ^1)- The taste breaking scales are related by 



2Ffm-^xi? 



(33) 



using the physical value for F in the last part. For lat- 
tice simulations with large volumes {V > 0.5 fm^) the 
e-regime observables will exhibit stronger taste breaking 
than the p-regime observables. This is the typical case for 
simulations where the dynamical quark masses are kept 
in the p-regime. For small volumes {V < 0.5 fm^) the 
relationship is reversed. However, in this case one might 
find that the dynamical mass is also in the e-regime so 
that the parametrization of Eq. ((27)) does not apply. Be- 
low we will examine the properties of the SRMT in both 
the weak and strong TB limits and explore the transition 
region between the two. 



VI. WEAK TASTE BREAKING 

A. Partition function 

Since one copy of the staggered Dirac matrix actually 
produces four tastes, it is common in simulations to take 
a fractional power of the quark determinant to produce 
the desired number of flavors in the continuum limit. We 
now consider the SRMT partition function with optional 
"rooting" given by 



Z'^f''^{{m},{n})= d[V]l[det[D 

fc=i 



TOfe 



|nfe/4 



(34) 



where the integration measure is over all the Gaussian 
weights and the powers can be either positive or neg- 
ative to produce a partially quenched theory. 
If we expand a determinant to order a^, we get 

det[(2?o + ruk) «) I4 + aT]""-^* ~ 

det[Vo + mkr''[l-a^{nk/8){Sl)] (35) 



with 



and we have used the fact that (Sk) = 0. The O(a^) 
partition function is then 

/ dm I n + m,r I 1 1 - E ^ ('5') I ■ (37) 

The last term in braces is the correction term for the 
partition function due to the TB. One can easily perform 
the Gaussian integrations over the taste breaking terms 
in the RMT to obtain the correction factor 

k=l 

+<3 {{Vo + mfc)-2) + ti ([75(2?o + mfe)]-') 1(38) 



with the dimensionless coefficients 
ti = a^V{4C3 + ACi 



t2 = aV(4C3 + 4C4 
t3 = a^V{C+ + C+) 



c\ 



ece) 

6^6 ) 



(39) 



Note that this expression can diverge as nik — )■ if 
i^Uk = 1. In general, this expression is not valid at very 
small masses since Sk can grow large, though for large 
enough masses it should be a good approximation. One 
could produce an alternate expression that is valid even 
at mil = using eigenvalue perturbation theory. Its con- 
struction will be outlined below. 

The correction factor can be further simplified using 
the identities 

(75(2^0 + ruky^) = v/mk 
([75(I?o + mfc)]-2) ^ ((Po-f mfe)-Vmfc) . (40) 

One still needs to integrate over the random matrix in 
I?o; however, all quantities can be calculated from known 
results for the partially quenched partition functions. 

If we write the partition function without taste break- 
ing in (|57|) as Ziq^[{m}, {n}), then the terms with a sin- 
gle trace can be readily evaluated as derivatives of the 
partition function without taste breaking using the sub- 
stitutions 



d,r.,ZNS{m),{n})^ J d[Vo]\D\nk{{Vo + mk)-') (41) 
dl^Zj,^{{m},{n})= (42) 
d[Vo]\D\ [nl {{Vo + 77ifc)-i)' - Uk {{Vo + mk)-^) ] 



where |-D| is short for the product of determinants in ([57)1 . 
This leaves only the two-trace term to be evaluated. This 
requires adding an extra quenched pair of quark species 
to evaluate the extra trace. 



(36) 



dm^dmjZN^+2{{{m}, ruf, rrib}, {{n}, 1, -1}) 



d[Vo]\D\nk{{Vo + mkry 



m f —niiy—rnk 



(43) 
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It is straightforward, though somewhat tedious, to 
evaluate these expressions in the microscopic limit. As 
an example, the one flavor partition function in the mi- 
croscopic limit is given by 



Zi(m) = /^(to) 



(44) 



while the partially quenched Nf — 2,Nb = I partition 
function can easily be written as a determinant of Bessel 
functions [2i,[23l, 



I^{m) m/,,+i(m) 7^2/^+2 (m) 
Ku{rhb) -mbKu+iirhb) mlK^+2{'mb) 



Z21 



f 



From these expressions we get that 

,2 \ / 



2i 



SRMT 



(45) 



l(to) + Iu+i{m) 



2 m 



+ (^1+^3)51-^3 



Iv-2{'m) + 21^{m) + I„+2{m) 



(46) 



where si is the complicated expression 
■'2 2iy 



-/^(m) + - 



'1/+1 



(to) + mK^+i{m)I^^^{m) 



+K4m) vilim) - {iv + 1)/'+i(to) 
2v{v + l) 



I,{m)I,+iim) (47) 



due to the two-trace term. 

From this, an expression for the quark mass depen- 
dence of the chiral condensate in the microscopic limit 
can be obtained by 



i;i(m)/So = drii InZf 



SRMT 



(rfi) 



(48) 



One can obtain expressions for any number of flavors 
through a similar procedure. These formulas would apply 
to lattice simulations performed in the e-regime, where 
TOjri ^ 1. For lattice simulations where this doesn't 
apply, one can instead consider observables as a function 
of a valence quark mass that is in the e-regime. This will 
be explored next with the quenched condensate. 



B. Quenched condensate 

One of the simplest valence observables one can look 
at is the quenched condensate. This can still be useful 
for comparisons with simulations of full QCD in the case 
that the valence quark mass is much smaller than the 
dynamical masses, so that the heavier masses will simply 
appear quenched compared to the light valence quark. 

The calculation simplifies considerably if we don't try 
to calculate the corrections to the quenched partition 



function first, but instead directly calculate the correc- 
tions to the quenched condensate from the definition 



So 



a.lni|«*^^({x,2/},{l,-l}) 



(49) 



Applying this to ([37)) gives 
2'x ti + Zt 



-'xxy 
,2 



2U 



xyyj 



2Zx 



(50) 



where Zx — dxZii{x,y)\y=x and similarly for higher 
derivatives. Using the expression for the microscopic con- 
tinuum partition function with Nf = Nb — 1, 



^11 = 



Iy(x) Xh+l{x) 

-K^{y) yK^+i{y) 



(51) 



we can now evaluate the derivatives. The result without 
taste breaking {Zx) is [l!] 

X [h{x)K^{x) + Iy+i{x)K^^i{x)] + v/x . (52) 

If we write the final expression with TB as Zx +X]fc ^fc-Zfc, 
then we have 

zi = Z2 = -2v'^/x^ (53) 
Z3 = 24 + 2/x - 2K,{x)[h+i{x) + I,-i{x)] (54) 

„3 



Z4 — ~2Ki,^i{x)I^+i{x)/x — 2vlx 



(55) 



Note that the single-trace TB terms enter only through 
the combination ^1+^2 oc C3 + C4, so that C\ and Cg 
don't contribute while C3 and C4 only contribute for v 7^ 
0. We will see in numerical simulations below that the 
quenched condensate is indeed fairly insensitive to small 
values of C4 for z/ = 0. 

The (partially) quenched condensate can also be used 
to calculate the eigenvalue density by inverting the 
Banks-Casher relation f29] 



E(m) 



1 
Y 



P(A) 

TO + iA 



■d\ 



The eigenvalue density is obtained from [Toj 

S(-a + e)-S(-iA-e) 



p(A) = lim 



27r 



(56) 



(57) 



Care must be taken when evaluating I](— lA— e) to use the 
correct formula for arguments with negative real parts. 
We can avoid this by making use of the fact that the 
condensate is an odd function of to (since p is even) to 
get 



P(A) 



i;(a) + i](-iA) 
2^ 



(58) 



Applying this to the quenched condensate in the micro- 
scopic limit gives the microscopic quenched eigenvalue 
density. The density without taste breaking is (30l | 



A 



[J,(A)2- J,+i(A)J,_i(A)] 



(59) 
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FIG. 1. Quenched chiral condensate for SRMT with A'' — 400 and v — for a range of values of C4 in the weak (left) and 
strong (right) taste breaking regimes. The condensate is fairly insensitive to C4 in the weak TB regime. In the strong TB 
regime it agrees well with the chiral GUE result (until finite effects set in at large x) with a scaled low energy constant Eq 
given in (f5^ . 



and the correction term due to taste breaking is 
hMX) [J.-i(A) - J.+i(A)] - {t, + U)^''^^^'^-+'^^^ 



X 



(60) 



Note that the ti and ^2 terms don't contribute at this 
level. This expansion will also break down close to zero, 
similar to the condensate. Expressions that are valid near 
zero can be obtained from eigenvalue perturbation theory 
as discussed in the next section. 



C. Eigenvalues 

The full staggered RMT is fairly complicated to work 
with; however, we can obtain relatively simple expres- 
sions for the effects of taste breaking on the splitting 
of the eigenvalues from eigenvalue perturbation theory. 
Since the spectrum without taste breaking contains de- 
generacies, we must take these into account. If is a 
set of degenerate eigenvectors for the eigenvalue iXk, 



(61) 



then the eigenvalues including the leading order pertur- 
bation are given by the eigenvalues of the matrix 



Ek=iXk + a-^lT'^k ■ 



(62) 



For simplicity, we first perform a unitary similarity trans- 
formation on Vq to put it in the form 




U^VoU= I iA 
iA 



(63) 



with A a positive diagonal matrix of the nonzero singu- 
lar values of W. The upper left block of zeros is of size 



h' X v, representing the zero modes, while the other two 
blocks along the diagonal are of size N x N. The above 
transformation can be absorbed into the taste breaking 
terms without changing their form, and we will not ex- 
plicitly write it anymore. The nonzero eigenvalues of 2?o 
written as above are simply plus or minus the diago- 
nal elements of A. In this basis the eigenvectors take the 
form 














In 


In 





In 


-In 



(64) 



We first examine the sector of the zero modes for 
which we define k — 0. The matrix Eq is just a projection 
onto the upper left i/ x i' block for all tastes of the taste 
breaking term. This gives the form 

«(^3ai + bvn + cy^) <S)S,f, + {Aif, + hAf, + CA^,) ® Ca^s (65) 

with implied summation over /i. The ^'s are v x v Her- 
mitian matrices and the 6's and c's are real scalars (as 
given in Appendix [B]). The terms are labeled accord- 
ing to their origin in the full SRMT and have the same 
Gaussian weights as the corresponding terms in the full 
SRMT. Here we see that only the pseudoscalar and ten- 
sor terms don't affect the splittings of the zero modes. 
Also, if we ignore the scalars and were to set C3 = C4, 
then the would-be zero modes are described by a chiral 
RMT, for which the eigenvalues are readily found. How- 
ever, typically we have that C4 ^ C3 so that, to a good 
approximation, we can set C3 — 0, which will give a dif- 
ferent distribution for the eigenvalues. 

We now examine the splitting within a degenerate 
quartet of eigenvalues. Here the 4x4 splitting matrix, 
^'iT^'fe, has the form 

idi + ddf,^ (g) + ixvp. ® + XAf, <S) S,f,5 (66) 
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FIG. 2. Quenched chiral condensate for SRMT with iV = 
400 and f — 1 for a range of values of C4 in the weak taste 
breaking regime. 



with = d3ij_+bvfj. and XAfi = d^fi+bAfi- All the d's are 
real scalars; again they are labeled according to their ori- 
gin, and their weight is exp{—aN^dx /VCx)- Since the 
^1,3,4,6 terms came from matrices in the original SRMT, 
they are different for each quartet. However, the bA,v 
terms came from scalars in the SRMT, so these variables 
are in fact the same in all quartets, and also in the zero 
mode sector. It is these terms that produce the leading 
order eigenvalue correlations between the different quar- 
tets and also allow the splittings within each quartet to 
be sensitive to the topological charge in the weak TB 
regime. Interestingly these are the terms in staggered 
chiral perturbation theory which don't contribute at one 
loop. 

In Sec. IVI Al we gave an expression for the partition 
function that could become singular as m — 0. Using 
the eigenvalues from perturbation theory we can derive 
a similar expression that treats the zero modes exactly. 
Since the 0{a) correction to the quark determinant van- 
ishes, we need to go to second order in eigenvalue per- 
turbation theory. The expression for the eigenvalues in 
second order degenerate eigenvalue perturbation theory 



VII. STRONG TASTE BREAKING 

Earlier we showed how the fermionic SRMT partition 
function maps onto the staggered chiral Lagrangian for 
weak taste breaking. There we took the large TV limit of 
the SRMT while scahng the terms a^a//? 1/N. This 
means that the a^VCx terms are kept fixed a,s N,V ^ 
00. Now we consider a different limit where the C'x are 
kept fixed so that a?a/j3 stays constant in the large N 
limit. Here the taste breaking must be included in the 
saddle point equations when going from the sigma model 
to the chiral Lagrangian, as mentioned in Appendix |^ 
The inclusion of TB breaks the Goldstone manifold from 
SU(4) down to U(l). Note that if only Ci 7^ then the 
symmetry is not fully broken down to U(l), but since that 
is not a likely scenario for common staggered actions, we 
will not consider it further. 

The result for the effective chiral Lagrangian in the 
strong TB regime is 



C 



1 



(68) 



where U € U{1). This has the form of a single fermionic 
flavor with a modified condensate given by 



S' = 



4So 



(69) 



with a'^D = a^VC/N. The factor of 4 is due to the four 
tastes, and the TB terms serve to decrease the effective 
condensate. 

Note that this result differs from what we would ob- 
tain from the large TB limit of the weak TB Lagrangian. 
That is, if we start with the Lagrangian (U) and take the 
saddle point in the non-Goldstone modes as — > 00 (or 
equivalently keeping the C's fixed as — >■ cxd), we would 
get the same form for the effective Lagrangian as in ([55)) 
but with a condensate that is simply Eg — 4So. This is 
because in the weak TB limit D would be scaled to zero, 
while it remains finite in the strong TB limit. Thus, in 
general, the two limits should be considered distinct even 
though some observables may exhibit a smooth transi- 
tion between them. We will explore this transition from 
weak to strong TB through numerical simulations of the 
SRMT. 



We can obtain an expression for the quark determinant 
(and hence the partition function) which is accurate to 
order by multiplying the determinants of the E^. If 
the determinants for each quartet were expanded in a, 
we would obtain an expression for the partition function 
identical to (1371) . However, to produce an expression that 
is valid at to = 0, the determinants must be handled 
exactly. We will not deal with that here. 



VIII. NUMERICAL RESULTS 

Here we present some numerical results obtained from 
averaging over 10,000 random matrix configurations ob- 
tained with N = 400 (and setting a = 1) at both v — 
and I' — 1. We varied the value of c = a^VC4 from to 
100, and all other TB coefficients were set to zero. 

In Fig. [T] we plot the quenched condensate from the 
SRMT. In the left panel, the taste breaking parameter 
c is varied between and 0.2. For c = we find the 
expected agreement with the plotted result from the chi- 
ral GUE (chGUE), Eq. for small x. For larger x 
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FIG. 3. Number variance for the SRMT with — 400 and 1/ — 0, left {v = 1, right), for a range of values of C4 showing the 
transition from weak to strong taste breaking. 



we see a discrepancy due to the finite N of the SRMT. 
As N is increased, the SRMT results will rise to match 
the chGUE curve. As c is increased, we see some small 
changes in the quenched condensate for small x, but not 
for larger x. This is in agreement with our expectations 
from the expansion ([55]) since we don't expect it to be 
valid near x = 0. At c w 0.2 we start to see changes in 
the condensate at larger x, signaling the end of the weak 
TB approximation. 

On the right panel of Fig. [TJ we show the quenched 
condensate as it leaves the weak TB regime and goes 
to the strong one. The scaled chGUE results shown are 
given by (normalized to a single flavor) 

tf°''e^x) = G™(sa;) (70) 

with s = Sq/So and using Eq. (IMl) . For c > 1 we see 
very good agreement between the SRMT results and the 
above formula, until the finite N effects set in at large x. 

In Fig. [5] we plot the quenched condensate v = 1 
in the weak TB regime. The chiral GUE result has an 
explicit V /x divergence which is matched by the 6*4 = 
data. As C4 is increased, the expansion of the conden- 
sate gives a TB correction proportional to —v^jx^ [Eq. 
([53|]. This corrected GUE result agrees with the data at 
c = 0.003 for X > 0.5, while a.t c — 0.01 the agreement is 
found only for approximately x > 2. For larger c the ex- 
pansion breaks down. At large c > 1 the 1^ = 1 quenched 
condensate enters the strong TB regime and will match 
the strong TB result of the 1/ — sector. 

As previously mentioned, the chiral condensate is re- 
lated to the eigenvalue density. Next we will look at the 
number variance which is related to the two-point eigen- 
value correlation function. It is defined as 

V^in) = {{[nin)-nf)) (71) 

where ((•)) denotes the ensemble average. The function 
n{n) is the number of eigenvalues between and £ in 



a particular configuration with £ chosen such that the 
interval has on average n eigenvalues. 

In Fig. [3] we plot the number variance for a range 
of values of C4, with all other coefficients set to zero at 
both ^ = and v = 1 topological charges. At C4 = 
the number variance matches that of the chiral GUE 
(T^*^^™) after the appropriate scaling 

yf-°(n) = 16K^'^°UE([- _ 2^]/4) (72) 

due to the fourfold degeneracy of the spectrum and the 
zero modes. Note that at C4 = there are really Aiy 
exact zero modes; however, we are considering this to be 
the limit of C4 — 0+, where for small C4 > the zero 
modes are split into 21/ positive and negative near-zero 
modes. This gives a shift of the number variance by 2i>. 

As C4 moves away from zero, we see the number vari- 
ance drop while retaining the same pattern of oscillations 
until a?VCn « 0.3. Here the large oscillations have essen- 
tially vanished. Upon increasing C4 the number variance 
continues to decrease and a set of smaller oscillations ap- 
pear that coincide with the result from the one flavor 
chiral GUE without any scaling. 

In Fig. |4] we can see the transition to the strong 
TB regime from the same data. As C4 is increased a 
larger range of the data falls on top of the chiral GUE 
curve. The point where the data starts to deviate from 
the chGUE curve is roughly equivalent to the value of 
a?VCi. For example for a?VCi — 10, the number vari- 
ance agrees with the chGUE result up to around n = 10, 
above which it begins to grow larger than the chGUE 
curve. This is similar to what one sees around the Thou- 
less energy [23j HE]; however, in this case it does not 
signal the breakdown of the RMT, but is instead the sig- 
nal of the restoration of an explicitly broken symmetry 
(taste). The behavior is independent of the value of N . 
If one observed a similar behavior of the number variance 
in lattice simulations, this could provide an independent 
estimate of the size of taste breaking. 
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FIG. 4. Number variance for SRMT with = 400 and u — 0, left [v = 1, right), for a range of values of Ca in the strong taste 
breaking regime. 



IX. SUMMARY AND FUTURE WORK 

We have presented a completed chiral random matrix 
theory for staggered fermions that includes all the taste 
breaking effects at order a? . The SRMT has been shown 
to be equivalent to the zero-momentum staggered chiral 
Lagrangian in the appropriate limit. We have also identi- 
fied a strong taste breaking limit where the SRMT maps 
onto a one flavor effective Lagrangian. 

The quenched condensate and number variance evalu- 
ated in the SRMT, clearly show the transition from the 
weak to the strong taste breaking regime. The number 
variance is particularly interesting since for strong TB, 
the TB scale a^VC seems to be directly related to the 
range of agreement to the one flavor chiral GUE result. 
We now need to test the results of the SRMT against 
lattice simulations with staggered fermions. 

One could also extend the SRMT to include an imagi- 
nary quark chemical potential. This can then be used to 
obtain a measurement of the low energy constant F j32r - 
Isg} and could potentially provide more information on 
the taste breaking. Similarly one could add a real chem- 
ical potential and study the effects of taste breaking on 
the complex eigenvalues [37|, |38|] . Additionally, one could 
explore the problems associated with rooting when the 
quark determinant is complex [39l . |40| . 



Appendix A: Mapping the SRMT to the staggered 
chiral Lagrangian 

The procedure for mapping the chiral RMT onto the 
zero-momentum chiral Lagrangian is by now standard 
0, mi- Here we will consider the case of only fermions, 
and will leave the addition of bosonic quarks (ghosts) for 
later. 

We introduce a set of Grassmann variables ijjf' and 



-0** , with t the taste index, s = ± the chiral index, and i 
the random matrix index. The determinant can then be 
written as 



det(I?o + m + aT)^ J d'ipd^e'f'^'^°+"'+'''^^'f' (Al) 



In this form the Gaussian integrals over the random ma- 
trices can be readily performed. This leads to a set of 
four fermion terms in the exponential. For the terms with 
scalars instead of matrices, we do not need to perform the 
Gaussian integrals now, and can wait until after expand- 
ing around the saddle point below. We will consider the 
different types of taste breaking terms separately, and in 
each case the Gaussian measure is taken to be 



exp 



{~I3N{T^T) /4) . 



(A2) 






iX \ 




a" 




J 




I3N 


( iA 










M 




\ 


iB ) 




I3N 



The four fermion terms generated by the two types of 
terms with random matrices considered can be summa- 
rized as follows: 



{ij'+ij^+Vkji^''-i,'-Tu) (A3) 



s=± 

These can be transformed into fermion bilinears by the 
Hubbard-Stratonovich transformation, yielding 

exp(-;3iV (pV) + a^j+pT'ip+ + a^)' p'^Tip-) (A5) 
exp(-;3iV {p\ + p^_) + a{l>+p+T%l}+ + aV^^p-FV'") (A6) 

respectively, where p is a 4 x 4 complex matrix and p± 
are 4x4 Hermitian. The Grassmann integrals can now 
be performed yielding a determinant. 
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In the absence of taste breaking, the RMT at this point 
would be 



da dct(a + 771)^+" det((7t + m)^e-"^('^^'") (A7) 



with (T a 4 X 4 complex matrix. In the large N limit, 
keeping ^/amN fixed (the microscopic limit) one finds 
the saddle point solution of cr = U/y/a, where U is a 
unitary matrix. The partition function expanded around 
this solution then becomes 



dU detiuy exp(^/^m7V ([/ + C/^)) . (A8) 



Matching to the chiral Lagrangian we get ^/a ~ 
T.0V/2N. 

To include the TB terms we must first choose how to 
scale those terms with N. If we had scaled the matrix p 
to have the same Gaussian weight as a, then we would 
get for ||X5|) 



exp 



{-aN {p'' p) + a^/^[tP+ pT'iJj+ + iIj- p'iT'iJj-]).{A9) 



The scaling of the term a a/ a/ f3 determines whether we 
are in the weak or strong TB regimes. If it is held con- 
stant then the term with p in the determinant will have 
a similar magnitude as the term with a and we must 
include it in the saddle point equations. 

For this strong TB regime, at the saddle point we find, 
for the term given in (jASp . 



t-i 



and in (IA6)) . 



p = (a//3)ra 



p+ = (a/2/3)ra-i 
p- = (a/2^)rcrt-i 



(AlO) 



(All) 
(A12) 



The two-trace TB terms do not contribute to the saddle 
point. The saddle point solution for a has the form 



a = ccxp(i^5 



(A13) 



for some c. The resulting effective Lagrangian for the 
strong TB limit of the SRMT is given in dBS]) with U = 
exp{i9). 

For weak taste breaking we can take a-y/a//? ^1/ y/N. 
Then the saddle point solution is unchanged by taste 



breaking and the TB terms can be expanded around the 
saddle point solution. Then the remaining Gaussian in- 
tegrals can easily be performed, yielding the set of terms 
given in Table HI 



Appendix B: SRMT 

For completeness, we will explicitly write down all the 
taste breaking terms from the full staggered RMT. The 
SRMT Dirac matrix is given by the form in Eq. (|13|) . 
The taste breaking terms T are given by the sum of the 
following terms; 




Ti 



'v 





The Gaussian weights for the individual terms are 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
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